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Solution of coupled vertex and propagator Dyson-Schwinger equations in the scalar 

Munczek-Nemirovsky model 

W. Detmold 

Department of Physics, University of Washington, Box 351560, Seattle, WA 98195, U.S.A. 

In a scalar ip 2 x model, we exactly solve the vertex and (p propagator Dyson-Schwinger equations 
under the assumption of a spatially constant (Munczek-Nemirovsky) propagator for the x field. 
Various truncation schemes are also considered. 
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I. INTRODUCTION 

The Dyson-Schwinger equations (DSEs) [l| are an in- 
finite system of coupled integral equations for the Green 
functions of a quantum field theory. In principle, their 
solution contains all possible information in the theory. 
However, the DSE for an n-point Green function involves 
various other Green functions of order m < n + 2. So to 
make progress, the system must be truncated at some 
stage and models must be made for the unknown Green 
functions. For example, in QCD, the DSE for the quark 
propagator involves not only the quark propagator, but 
also the gluon propagator and the full quark-gluon ver- 
tex, both of which need to be calculated or modelled in 
some way for the equation to provide a solution for the 
quark propagator. 

Notwithstanding the need for model input, the DSEs 
provide a relativistically covariant, non-perturbative ap- 
proach to field theory that has enjoyed considerable suc- 
cess in recent years. In QCD, models based on the DSEs 
are able to describe accurately many facets of meson 
and, more recently, baryon physics in free space and at 
nonzero temperature and density 0, Additionally, 
recent DSE studies are now providing insight into the 
infrared behaviour of QCD Green functions and the na- 
ture of confinement 0,13- Sophisticated studies have also 
been made of the DSEs in QED and other theories @>E1- 

Almost all DSE based studies to date focus on the 
equations for some of the two point functions of the 
theory, making various models for the remaining two- 
and three-point functions. In this paper, we examine 
a simple theory and solve a larger system of equations. 
We consider a ip 2 x scalar field theory with a Munczek- 
Nemirovsky (MN) model propagator for the x particle 
(which we introduce in Sec. [Hj) and solve the coupled 
DSEs for the tp propagator and 3-point vertex. An ex- 
act solution is possible in this model because relations 
amongst the Green functions allow the remaining system 
of DSEs to close without truncation. We also investigate 
the effect that imposing various truncations has on the 
solutions. 

The central result of this analysis is that there is no 
real solution to the untruncated model unless the cou- 
pling is imaginary despite the model having solutions for 
real couplings when additional truncations of the DSEs 
are made. This is shown by the divergence of the dia- 
grammatic expansion of the vertex. 

The rest of this paper is structured in the following 



manner. We introduce the details of the model in Sec. ITTI 
and the various truncations we employ in Sec. IIIII In Sec. 
IIVI we show that the DSEs of the Munczek-Nemirovsky 
model close and derive their exact solution where it ex- 
ists. Finally, we present our conclusions in Sec. [V] 



II. SCALAR MODEL 

We consider the theory of two scalar fields, <p and \i 
described by the interaction Lagrangian: 



Ant = 9 V 2 X • 



(1) 



and work in Euclidean space. This scalar theory is pri- 
marily studied as a precursor for an investigation in 
QCD, but is also interesting in its own right and is some- 
times used as a simplified description of meson-nucleon 
interactions. We define the momentum- (and position-) 
space propagators for the fields ip and x as S(p) (S(x, y)) 
and D(p) (D(x, y)) respectively. These objects satisfy the 
Dyson-Schwinger equations: 



and 



D- 1 (p)^D^ 1 (p) + n(p). 



(2) 



(3) 



where the bare propagators are Sq 1 (p) — (p 2 + m 2 ) for 
the ip field of mass m and Dq 1 = (p 2 + M 2 ) for the x 
field (mass M). The respective self-energies are defined 



as 



and 



£(p) = J g D(q)S(p-q)f(p-q,p) 



(4) 



^{p) = \J^9 S(q)S{p + q)f(q,p + q), (5) 

where T(p, k) is the proper three-point vertex. (Here and 
in what follows momentum integrals are abbreviated as 
J = J (2tt) 4 •) These equations are shown diagrammati- 
cally in Fig. An additional (constant) tadpole contri- 
bution to Eq. (@J is omitted as it can be absorbed into the 
definition of the renormalised mass, m 0- Field renor- 
malisation constants are suppressed as the ultraviolet di- 
vergences of the theory will be eliminated by the specific 
form of the x propagator that we introduce below. 
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FIG. 1: The Dyson-Schwinger equations for the (a) tp and 
(b) x propagators and (c) vertex function of the scalar theory 
of Eq. 0. Solid (dashed) lines are tp (x) propagators and 
grey circles denote full one particle irreducible components. 



the x field is constant throughout space-time. In momen- 
tum space this results in 

g 2 D(p) = (2^) 4 a S\p) , (7) 

where a is the (constant) strength of the correlation. 
This is analogous to the approximation of Munczek and 
Nemirovsky (MI\ Q I12| which has often been applied to 
model QCD [H Il4l It was originally motivated 

as a simple means with which to provide the necessary 
infrared enhancement in the kernel of the quark DSE of 
QCD. However its physical foundations are fairly tenuous 
as it does not allow momentum transfer between quarks 
or their consequent scattering. Whilst it has this and 
other difficulties 0], it also has the distinct advantage 
that the various integral equations of the DS system re- 
duce to algebraic equations and it remains an interesting 
model. 

With this scalar MN model in place, we will explore 
varying possible truncations in the following section and 
then show that they are unnecessary by solving the 
full set of equations and calculating the exact solutions 
(where they exist) in Sec. IIVI 



The dependence of these equations on the vertex func- 
tion r means that they do not form a closed, computable 
system. An often used approach is to truncate the sys- 
tem at this level by introdu cing a model for T and solve 
the resulting equations [sHj^llOf. 1 We investigate various 
possible truncations in Sec. IIIII 

However, the vertex function satisfies its own DS equa- 
tion [an inhomogeneous Bethe-Salpeter equation — see 
Fig. Wi c )] m terms of the propagators and an additional 
unknown function, the fouv-tp scattering kernel, K: 

r(p,fc) = g^t&k) (6) 
= 1+ [ K(p,k,q)S(p + q)T(p + q,k + q)S(k + q). 

Jq 

Here we have rescaled T — ► T to absorb the bare cou- 
pling of the theory. The scattering kernel, K, occurring 
in this equation is two particle irreducible with respect 
to the tp legs and contains no annihilation contributions 
0. The same kernel also enters in investigations of two- 
body bound states in the Bethe-Salpeter approach 
By keeping Eq. |JBJ, it is possible to truncate the system 
at a higher level, modelling K rather than T. This auto- 
matically ensures that the Dyson-Schwinger and bound 
state equations are truncated consistently. 

In this paper, we will assume that the propagator for 



In phenomenological approaches to QCD, only the analogous 
equation for the quark propagator is retained and the com- 
bination of the gluon propagator and quark gluon vertex is 
modelled Ij. 



III. TRUNCATION SCHEMES 

In this section, we shall consider a number of trunca- 
tions that can be applied to solve the DSEs above, each 
of which successively improves on the last. For simplic- 
ity, we omit the momentum dependence of propagators 
and vertices in what follows since the MN model decou- 
ples all differing momenta. In fact, the momentum only 
enters into the propagator and vertex DSEs though the 
bare propagator, Sq 1 — p 2 + m 2 = x, so T, S, S and K 
depend only on x. 

It is possible to define the dimensionless quantities 
S' = y/\a\S, £' = \/MX r ' = r, x' = x/y^\ and 
K' = K/\a\. If one works in terms of these new objects, 
all dependence on the magnitude of a will factorise in our 
analysis. Consequently, only the sign of this parameter 
is important. In order to avoid a proliferation of factors 
of a/ 1 a |, we shall use the original (unprimed) quantities 
and only consider a — ±1. 

A. Rainbow truncation 

For this scalar theory, the rainbow truncation is defined 
by 

r r ainbow(P, k) = 1 . (8) 

In QCD, when the analogous truncation of the quark 
DSE is combined with the ladder truncation of mesonic 
Bethe-Salpeter equations, the resulting models success- 
fully describe the pseudo-scalar and light vector me- 
son spectrum 17] with additional higher-order correc- 
tions providing only minor modifications 
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B. One- and two-loop truncations 

In the algebraic MN model, it becomes simple to in- 
clude any finite order of contributions to the vertex. As 
no momentum flows along \ propagators and there are no 
complicating Dirac structures, every diagram at a given 
order contributes the same factor, For example, the one- 
and two- loop vertices are given by 



Ti loop = 1 + a S 2 



r 



2 loop 



l + a S z + Aa 2 S i , (10) 



respectively, where the factor of 4 arises from diagram 
counting (see Sec. II V A")) . These additional contributions 
to the vertex are not perturbative corrections; they in- 
volve the full propagator functions. Inserting these forms 
into the propagator DSE leads to 4th- and 6th- order 
equations for S whose physical solutions are also shown in 
Fig. [3 as the short- and long- dashed lines, respectively. 
The lower panels of the figure show the corresponding 
self energies and vertex functions. For comparison, we 
also show the solution (dotted line) when only the ladder 
dressings of the vertex are included to second order (that 
is, T2 loop with a factor of 1 instead of 4) which we label 

aS ^2 loop- 



Vertex ladder truncation 



FIG. 2: Dependence of the propagator (upper panel), self- 
energy (middle panel) and vertex function (lower panel) on 
the variable x = p 2 + m 2 for a = 1. Solutions are shown 
for the bare- (dot-dashed), the one- loop- (short-dashed), two- 
loop- (long-dashed), two- loop-ladder- (dotted) vertices and 
for the ladder truncation of the vertex DS equation (solid). 



the rainbow truncation is combined with the Munczek- 
Nemirovsky propagator, Eq. J7J), the DSE for the if prop- 
agator reduces to a simple quadratic equation whose so- 
lutions are 



If (instead of explicitly approximating the vertex) 
we make a ladder truncation of the vertex kernel, 
K(p,k,q) = g 2 D(q), we are selecting an infinite class 
of planar diagrams (with non-perturbative propagators) 
to include in the vertex. At each order, we include the 
diagram that recursively dressed the contribution with 
one less loop |l5|, [lSf. However, even at the two loop 
level, contributions are already being omitted (compare 
T 2 loop and T' 2 loop ). 

With this truncation, Eq. 10 reduces to 

T{p, k) = 1 + g 2 f D(q)S(p + q)T(p + q,k + q)S(k + q) . 

Jq 

(11) 

Using the MN propagator, Eq. J7J), gives the solution 



S(x) 



x ± \/x? 



+ 4a 



2a 



(9) 



laddc 



00 1 

Y> S 2 Y = -^- 

^ K ' 1 - a 

i=0 



S 2 



(12) 



The physical solution 2 and the corresponding self energy 



are shown as a function of x 
curve in Fig. for a = 1 . 



p 2 + m as the dot-dashed 



provided a S 2 < 1. When this is inserted in the propa- 
gator DSE we arrive at a cubic equation for S: 



s- 1 = 



a S 



I- a S 2 



(13) 



2 It is important to note that for each truncation considered here, 
the algebraic equations have only one solution that asymptotes 
to the bare propagator for large momentum or particle mass. 



The physical solution to this equation is given for a = 1 
as the solid line in Fig. [3 The sum in Eq. (|12f> is seen to 
converge for all space-like and time-like x. 
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IV. COUPLED VERTEX AND PROPAGATOR 
EQUATIONS 



From comparison of the rainbow, one-loop, two-loop 
ladder (T' 2 loop ), and full ladder solutions of the previous 
section (see Fig. 01 , it is apparent that the approach to 
this ladder sum is monotonic and inclusion of such ladder 
diagrams beyond the first loop only results in small mod- 
ifications |l5j . However, further comparison with the full 
two- loop vertex solution (corresponding to T2 i 00 p) indi- 
cates that the omitted diagrams may be important. In 
the simple MN model defined by Eq. Q, it is possible to 
go further and consistently include all contributions to 
the vertex. This is the main result of this article. 

The generating functional of the theory we consider is 



z[J,j] 



(14) 



VtpVx exp 



£ - (f, J) - (x,j) 



where J(j) is an external source for the f{x) field, 
C is the sum of the interaction Lagrangian, Eq. JQl, 
and the usual kinetic pieces and the notation (a, b) = 
J dxa(x)b(x). From this we can make a Legendre trans- 
form and define the effective action, A[ip, x] (with if — 
—6W/5J and x — — SW/Sj), and thereby the one parti- 
cle irreducible vertices, 



p(n|m) 

*■ xi,...,Xn,yi, — ,yn 



sn+m 



(15) 



Stp Xl . ..Stp Xn Sx yi ■ ■■&X Vm 
The DSEs arise from appropriate functional variation of 



this effective action. As usual, we can identify the posi- 
tion space propagators with the inverse two-point, onc- 
particle irreducible (1-PI) vertices 0. For example, 



S(x,y;J) 



S 2 W 

SJ X 5jy 



S 2 A 



(rg ))" 1 



Thus, 



(16) 



K(w,x,y,z) 



A 2 r (2|0) 

r (4|o) _ 1 yz 

wxyz r c 

d(p w dip x 
5Y>{y,z;J) 
SS(w, x; J) ' 



(17) 



evaluated in the presence of external sources. For the last 
equality, we have used the position space analog of the 
DSE for the <p propagator [Eq. (J2J] to express S~ 1 (x, y; J) 
in terms of y; J), noting that SSq/SS vanishes. 

By momentum conservation, the full 1-PI four-point 
scattering kernel, K{p 1 q,r), is a function of three ex- 
ternal momenta. However, in the MN reduction of the 
theory, no momentum is transferred through x propaga- 
tors and the scattering kernel describes two separate mo- 
mentum flows. Consequently we only require K(p,q,0). 
Further, since we only probe the q = limit of the ver- 
tex function in the vertex DSE, the (p propagators in 
Eq. jnjl must carry the same momentum and we only 
require K(p,q — p, 0) = S~E(p)/SS(p). 

Inserting the self energy from Eq. (J2J and applying the 
derivative gives 



K(p,p,0) = 



SS( P ) 



SD(q) 

6S( P ) 



S(p - q)T(p -q,p)+ D(q)5\p - q - p)Y{p - q,p) + D(q)S(p - q) 



5T(p-q,p) 
6S(p) 



(18) 

In principle, Eq. J2J determines the functional dependence of D(q) on S(p) through the polarisation loop in Eq. J5J; 
pre- and post- multiplying both sides of Eq. © by factors of D(q) and taking the derivative gives 



= 2D(q)S(q - p)T(q - p,p)D(q) + D(q) J S(k) S -^-^S(q k)D(p) 



(19) 



However, the MN approximation, Eq. Q, presumably 
already invalidates Eq. so this result may be mislead- 
ing. As a simple alternative to Eq. Ijl9(l , we also consider 
SD(p)/SS{k) = 0. 

We proceed by applying Eq. Q to Eqs. JEJ and 
(JTSJ to give 



ST. 



r = 1 + KS 2 T . 



(20) 



(21) 



K = a 



r+ 4 



■K D . 



(22) 



If Eq. I|19|l is used to determine 6D/SS, there is an addi- 
tional contribution, 



K 



D 



l s r 



9 <sr 

2ST + S 2 — 
ob 



(23) 



to the kernel. For SD/5S = 0, Kd = 0. Momentum 
dependence is again suppressed in these equations as all 
quantities are at the same momentum. Combining the 
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equations for the vertex and the kernel gives 

r = 1 + a S 2 T (r [1 - 2a S 2 T] + Sjg [l - a S 2 T] 

for 5D/8S determined from Eq. (|19|l . and 

.ST 



(24) 



r = 1 + a S 2 T 2 + a S 3 T 



SS 



(25) 



for SD/SS = 0. These functional differential equa- 
tions are subject to the physical boundary condition 
that r[0] = 1. This is simply the requirement that 
for an infinitely heavy tp particle (whose propagator 
S = 1/m 2 "^-^5° 0), the vertex reduces to the bare vertex 
as loop corrections involving the if propagator vanish. 
Either of these equations, together with their boundary 
condition 3 and Eqs. J3J) and (|20|l close provided the func- 
tional derivative ST/SS can be determined. 4 



A. Diagram counting 

Physically, the use of the Munczek-Nemirovsky-like 
propagator, Eq. Q), means that x field loops dressing 
the propagation of the <p field or dressing the vertex carry 
no momentum. Consequently, objects such as the ip self 
energy and the vertex can be expressed as appropriately 
weighted sums of products of dressed tp propagators at 
the same momentum. That is, T and £ are analytic func- 
tionals of S and can be written as 



r = g/,(a s 2 y. 



i=0 



and 



S = a S Bi(a S 



2 V 



(26) 



(27) 



i=0 



As a corollary, ST/5S exists and can be defined. Indeed 
Eqs. i|24|) and 1251) can both be rewritten as recursions 
amongst the coefficients /j. For 6D/5S = 0, the relation 
is 



Si = S i0 + »•/ + l)/j/<-i-i 

3=0 



(28) 



These coefficients simply count the number of TPI con- 
tributions to the vertex dressing with i loops in terms of 



3 At x = 0, the propagator DSE reduces to T[S] = l/(aS 2 ) which, 
when inserted in the vertex equation [Eq. 1251 1 . implies S = 
i-y/i/a and rfS 1 ] = 77. However, here the functional dependence 
of the vertex is incompatible with the boundary condition and 
the solution is unphysical. 

4 A related method of functional closure wa s presented very re- 
cently in a perturbative context in Ref . |2(1 . 
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FIG. 3: Dependence of the propagator (upper panel), self- 
energy (middle panel) and vertex function (lower panel) on 



the variable x — p + m for a 



-1. Solutions are shown 



for the bare- (long-dashed), the one- loop- (dot-dashed), two- 
loop- (dot-dot-dashed) vertices, for the ladder truncation of 
the vertex DS equation (dotted), for the full solution with 
SD/SS — (short-dashed) and for the full solution with 
SD/SS from Eq. I|19^ (solid). Fo clarify the various solutions 
for the propagator, the shaded region is expanded in the inset. 



full propagator lines. 5 The first few coefficients are: 1, 
1, 4, 27, 248, 2830, ... and their asymptotic behaviour 
is f n ~ (2n + l)!!/e from which it is apparent that 
the sum in Eq. (|26|l does not converge (nor is it Borcl 
summable) for a S 2 > 0. For the case of 5D/5S given 
by Eq. (|19fl . an additional term 

i-2 i—j—2 

- 2 E E (*+i)/iA/i-i-fc-2 

j=0 fe=0 



5 Th e enumeration of this series is considered in Refs. L'lL l22l 12^ . 
I'24H . but no closed form for the i th element exists to the author's 
knowledge and the recursion of Eq. 1281 provides an efficient 
calculational method. 



6 



appears on the right-hand side of Eq. I|28|l and the series 
is similarly divergent. The self energy coefficients, e^, 
can also be calculated and a generating function can be 
constructed for them pEj . 

The apparent lack of solution is supporte d by the fail- 
ure of various standard numerical methods |2q to find a 
solution to either of the differential equations for F with 
a > 0. This result is perhaps not unexpected since the 
theories defined by the Lagrangian of Eq. and the re- 
lated ip 3 theory are unstable as they have no ground state 
H USE! EllHl However, the effect that the use of the 
MN model for the x propagator has on these conclusions 
is difficult to assess. 



Numerical solution for a < 



If we consider a < (as discussed at the beginning of 
Sec. 11111 the magnitude of a is unimportant), solutions 
of Eqs. (|24l )and exist. For this to be the case, the 
theory either has an imaginary coupling, g — ih (h G-ffiJ), 
or the correlation between the x fields at two points is 
negative so the fields are imaginary valued. The first sce- 
nario corresponds to a non-Hermitian Hamiltonian that 
is T-'T-symmetric. Such theories retain a positive energy 
spectrum and have been considered in detail by Bender 
et al. HU. 

The numerical solution of the vertex differential equa- 
tions and their boundary condition is straightforward for 
a < (the functional nature is irrelevant because of the 
decoupling of propagators at different momenta). Once 
we have determined the functional dependence of the ver- 
tex on S, we can combine that with Eqs. (|2J) and (|20|l and 
obtain exact solutions to the MN model. 

Figure depicts the momentum dependence of the ip 
propagator, the self-energy and the vertex as a function 
of x for a = — 1 and for the two alternate determinations 
of 8D/SS. The solid line corresponds to SD/5S deter- 
mined from Eq. l|19fl . while the short-dashed line is for 
5D/SS = 0. For comparison, we also show the (a = —1) 
solutions obtained using various of the truncations of the 
preceding section. For a < 0, most of the algebraic solu- 
tions become complex at some critical value of x, hence 
the abrupt termination of the solutions in the figure. Nu- 
merically we find that the full solution for 5D/SS = 
has a critical point x cr ;t = 1.518-y/|a|, although it is not 
known whether this critical behaviour has any physical 
significance. For SD/SS determined from Eq. (|19|l . the 
solution exists for all x. It is also evident from the figure 
that the various analytic solutions beyond the bare vertex 
provide a reasonable approximation to the full solution 
for moderate to large momenta as should be expected. 



C. Higher-point vertices 

We now know the functional dependence of T[S] (— 
r( 2 W) and, through Eq. 0, of K [S] (= T^ )) for 



FIG. 4: Behaviour of the higher point Green functions as 
functionals of S and functions of x as described in the text. 



a < 0. Consequently, we can (numerically) construct 
every 1-PI Green function of the theory with zero or one 
external x leg carrying zero momentum and any number 
of pairs of tp legs carrying the same momentum. Figure^] 
shows r( 2 l°) (short-dashed), T^ 1 *) (dotted), r( 4 l ) (long- 
dashed), r< 4 W (dot-dashed), r( 6 l°) (solid) as functionals 
of S (upper panel) and as functions of x (lower panel) 
given the numerical solution for S, all with a = — 1 and 
SD/5S obtained from Eq. (P|) . 

Since we do not know the behaviour of K(p, q ^ p, 0), 
we cannot apply the results of this analysis to the study 
of massive bound states within the model as was done in 
Ref. for the ladder vertex truncation. 



V. SUMMARY 

The analysis in this paper demonstrates that the ex- 
act MN model defined by Eq. has no solutions for 
positive values of the coupling a. However when addi- 
tional truncations are introduced for the vertex or for the 
scattering kernel, solutions can be found. (It is of course 
perfectly reasonable to define the MN model to include 
the bare vertex truncation.) For the case of a < we 
solve the full MN model and investigate how well various 
truncations approximate the full solution. 

The extension to a more QCD-like theory (which was 
the motivation for this study) through the use of fermion 
and vector boson fields appears possible after the compli- 
cations introduced by the Dirac structure are addressed. 
This will result in coupled, partial differential equations 
for the propagator dependence of the various Dirac com- 
ponents of the fermion-gauge boson vertex. Since QCD 
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and more realistic models of it do not suffer from the 
fundamental difficulties of the scalar theory considered 
here, the problems of solution existence that we encoun- 
tered would likely not persist. As a further extension of 
this work, there may be some scope to extend this anal- 
ysis to models with less restrictive \ propagators such as 
Gaussian forms [32| . 
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